Annals of the Academy of Romanian Scientists 
Series on Physics and Chemistry Sciences 
ISSN 2559 — 1061 Volume 5, Number 1/2020 


70 


ALPHA-RESONANCES IN MEDIUM 
AND HEAVY NUCLEI 


II. Review of microscopic models 


D.S. Delion*, V.V. Baran’, and A. Dumitrescut 


Abstract 


We review microscopic approaches describing a-like resonances in nuclei 
with Z > 50. We give a microscopic explanation of a-like rotational bands in 
the 4°Ca-+-a system in terms of single particle Gamow resonances. On the other 
hand, we show that for a-decaying nuclei the decay width can be described only 
in terms of a preformed a-cluster, existing on the nuclear surface in addition 
to the standard mean field cluster. Thus, we use a semi-microscopic hybrid 
model combining the mean field formation with a preformed a-cluster in order 
to explain the order of magnitude of the experimental decay width. Finally, we 
analyze proton-neutron versus a-like correlations. 
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1 Introduction 


a-particles were first evidenced by the discovery of the a-decay mode more than 
100 years ago, marking the birth of nuclear science. The a-decay mode occurs for 
medium/heavy nuclei with Z > 50. Its explanation was proposed soon after the 
establishments of quantum mechanics [1, 2] in terms of the quantum penetration of 
a preformed a-particle through the Coulomb barrier. The microscopic description of 
the a-particle preformation [3, 4, 5, 6, 7] still remains a theoretical challenge, due to 
the fact that a quartet structure, built in terms of a collective superposition between 
two-proton plus two-neutron single particle orbitals generated by a standard nuclear 
mean field, underestimates the experimental decay width by two orders of magnitude 
[8]. 

First of all, we will give an explaination of a-like rotational bands in the 4°Ca+a 
system in terms of the microscopic a-particle formation amplitude by using single 
particle Gamow resonances provided by the standard nuclear mean field. 

On the other hand, the analysis of decay data between ground states shows 
that the probability of the a-particle formation on the nuclear surface, estimated 
by dividing the decay width to the penetration probability, is proportional to the 
fragmentation potential, given by the difference between the Coulomb barrier and 
Q-value. This quantity does not enter the standard estimation of the preforma- 
tion amplitude in terms of single particle orbitals generated by some nuclear mean 
field. A simple two-body quantum calculation of Ref. [9] clearly shows that this 
proportionality can be explained by considering a preformed cluster in a pocket-like 
potential centered on the nuclear surface, i.e. similar to the a-like resonant state 
seen in light nuclei by scattering experiments. 

The aim of this paper is to review the microscopic theoretical approaches con- 
cerning the resonant character of the a-decay phenomenon. 


2 Microscopic a-formation amplitude 


In a microscopic description of the a-decay process the internal part of the a-core 
potential is replaced by a wave function, seen as a formation amplitude given by the 
following overlap between initial and final configurations [4, 5] 


F(R) = (aD|P) = | eodxp [Wa(Xa)Vp(xp)|* Up(xp) = fi) . (2.1) 
where by x we denoted the internal coordinates of the fragments. Therefore, the 
final result depends upon the relative radius between emitted fragments R. The 
antisymmetrisation a-core effects are neglected here, because we will estimate this 
quantity beyond the geometrical touching point (i.e. in the barrier region), where the 
two fragment have a small geometrical overlap and therefore Pauli effects are strongly 
diminished. The wave function of the a-particle wq is expressed as a Gaussian in 
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relative and center of mass coodinates, which is actually a product of the lowest 
harmonic oscillator (ho) states with the parameter [4] 


Ba = aren ~O.5fm- . (2.2) 


Theoretical nuclear structure is built on the single particle (sp) mean field concept. 
Thus, the wave functions of the parent and daughter nuclei in Eq. (2.3) are built 
as various superpositions of antisymmetrised products between sp orbitals (Slater 
determinants) generated as eigenstates of some nuclear mean field. The mean field 
is generated by using the standard Hartree-Fock (HF) selfconsistent procedure for 
many-body systems above magic nuclei, involving only few valence sp orbitals, and 
therefore the wave functions of even-even parent/daughter nuclei are Slater determi- 
nants built on these orbitals. For superfluid systems between magic nuclei one uses 
the Hartree-Fock-Bogolyubov quasiparticle method. A more pragmatic approach 
uses a parametrised mean field with a Woods-Saxon shape. 

In the absence of the a-core antisymmetrisation the wave function of the parent 
nucleus can be approximately factorized as a product between a four-body quartet 
state and the daughter wave function, i.e. 


Up =UyVp, (2.3) 


and therefore the formation amplitude becomes 


R) = | tts %a)Vo = (Hala) . (2.4) 


Typical examples of nuclei with an a-like structure are first of all double-magic plus 
quartet structures like *9Ca+na, !°Sn+na, 22 Pb+na. For n = 1, the procedure 
to estimate the overlap integral (2.4) consists in expanding the quartet in terms of 
products between pp, nn and pn pairs of sp orbitals. Due to the fact that the above 
mentioned systems are spherical we can use the angular momentum coupling scheme 


|a) = 2. 10) ’ (2.5) 


in terms of the quartet operator, given by the following superposition 


t ! T i) 
Qi,M, = da (Suu) ) [P},(@p) ® P}y(nm)] 
by 
+ FY AC) [PO @PLOM],,, + 28 


where J, is the quartet spiNgigenvalue index. The collective pair operators 


Per) = S>Xy, (7557 a, (7557'9) , (2.7) 


jy! 
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depend on the normalized basis pairs, built on sp orbitals 


a ® cl, 2 A (2.8) 


= sighah 
Crm 3 ) I) IuMy 


1+ Orgyrj" | 


where J,, is the pair spingigenvalue index. The isospin index is T = p,n and the sp 
orbitals ae generated by the nuclear mean field, depend on spherical quantum 
numbers j = (€,j”) (eigenvalue, spin and parity). The X-two-body coefficients are 
given by the equation of motion procedure with some two-body interaction V. 


2, Phy GaP +V, Pl yr) = By, (er Ph ag, (Cr). (2.9) 
By supposing that the basis pairs operators satisfy a boson commutation rule 
[Canga (1515 711) Ch, any (T2J257343)) = 51,2 ; (2.10) 


one obtains the standard system of the Tamm-Dankoff Approximation (TDA) 


1 . NS : . 
5 S oe 9 |V irk: a 8 OG A al Th) = [Zy,(77’) a ae Epigt | AY, (r9;7'7') , 
kk’ 
(2.11) 


for (tr') = (pp), (nn) and (pn) systems. The four-body Z-coefficients are also 
provided by the equation of motion procedure [10] 


Og | CLO ie (2.12) 


By using the equations for pairs (2.9), one obtains the following system for Z 


Ay, Ay Za) _¢, (Nu Ne Ze (2.13) 
Ay, Ao Zz? —le\ Nog Nos A 
where the norm and interaction matrices are defined as follows 
y 
= t top t on 
Nim = (Ol[PHa)ePl@K)] [Plime Pim |0) 
1 
Hem = 5 [E(k) + E(k’) + E(m) + E(m’)] Nem ; (2.14) 


in terms on spin-isospin pair index: 1 = (J, pp; Jinn) , 2 = (Jupn; Jpn). Their 
detailed form is given in Ref. [10]. 

The overlap integral (2.4) can be now estimated for each eigenstate [, (2.5) pro- 
vided by the system (2.13) in a standard way. Namely, the product between two 
proton and two neutron sp orbitals, entering the expansion (2.6), can be rewrit- 
ten in terms of relative and center of mass (cm) Moshinsky coordinates in the ho 
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Figure 1: The hindrance factors (2.17) as functions of the excitation energy for 
I} =07,27,47,67 four-particle eigenstates in the system *4Ti =*9 Ca + a [11]. 


yop “py? 


basis. Finally, the overlap between the a-particle wave function and the resulting 
ho functions in pp, nn and pn relative coordinates can be written as the following 
superposition 


F1,(R) = (bal Qt, |0) = S> Wr, (NoLa)® Yet, (R) - (2.15) 
NolLe 


in terms of radial ho functions depending on cm coordinate of the a-core system 
and four times the sp ho parameter 
_ Mywn 


p=——, fun = dia’ (2.16) 


Details are given in Section 5.1, according to Refs. [10, 12]. The nuclear structure 
information is carried out by the pair of X-coefficients (2.7) and quartet Z-coefficients 
(2.6). Let us mention here that, according to Eq. (2.15), the a-particle in some state 
I, is born from the ground state by the action of the quartet operator Ql. The 
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Figure 2: The formation amplitude for the system Pb+ a versus radius 


[10]. 


lowest eigenvalue J, = OF corresponds to a transition between ground states. The 
higher excited eigenstates are called a-like resonant states. 

These states can be probed in light systems by the elastic scattering of a- 
particles. We estimated their structure for “Ti=*?Ca-+a in Ref [11]. As sp orbitals 
were considered all bound and narrow Gamow sp resonances in “°C'a, generated 
by a Woods-Saxon potential with universal parametrisation. Then pp, nn and pn 
X-coefficients of the pair excitations were estimated, by using the surface-delta in- 
teraction in the TDA system of equations for 4?Ti, 42Ca and Sc, respectively. 
The doubly magic nucleus 4?Ca was considered as a vacuum state. Finally, the 
Z-coefficients and the a-formation operlap integral to the quartet eigenstates were 
computed. They are given in Ref. [11]. In Fig. 1 of the same reference we estimated 
the “hindrance factors” defined by the ratio 


HR = 2 217 
= Sit ? ( . ) 
beween spectroscopic factors 
Sr, =| \Fr,(R)/?RdR , (2.18) 
0 


to some excited state [7 and ground state OF. The positions of maximal values 
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R (fm) 
Figure 3: The equivalent local potential of the system 7/2 Po =2°8 Pb-+ a (solid line) 


and Coulomb interaction (dotted line) [10] versus radius. By a horizontal line is 
given the Q-value of the a-decay from 7!" Po [10]. 


revealed a very good description of the measured quasimolecular resonances, as can 
be seen in the Table 1. 
Table 1. Energies of a-like states in “+ Ti [11]. 


Eexp (MeV) | Eth (MeV) 


The microscopic structure of a-resonances was also estimated for 7!2Po=?°8Pb+a 
in Ref. [10], were this microscopic method was initally proposed. The low-lying 
states in 7!°Po, ?!°Pb and ?!°Bi were described as collective superpositions of pairs 
by using the surface-delta interaction in the TDA system of equations. The double 
magic nucleus ?°°Po was considered as a vacuum. We considered sp bound states 
and narrow Gamow resonances, generated by the Woods-Saxon mean field with 
universal parameterisation. Theoretical predictions for energies are in a reasonable 
agreement with experimental data, as can be seen in Table II of Ref. [10]. In spite 
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of its simplicity, the agreement of energies for low-lying states in ?!*Po was good, as 
can be seen in the Table 2. 


Table 2. Energies of low-lying states in 7? Po [10]. 


Eexp (MeV) | Ein (MeV) 


Several high-lying a-like states with relative large hindrance factors were pre- 
dicted to be seen in future experiments. 

The main advantage of this heavy system is the existence of the a-decay phe- 
nomenon. The formation amplitude for transitions between ground states I, = 07 
is given in Fig. 2. One notices the peak on the nuclear surface denoting that the a- 
particle is indeed created at small densities. By using the formation amplitude as a 
wave function satisfying some Schrédinger equation, is is possible to estimate the lo- 
cal equivalent potential corresponding to this solution. It is given in Fig. 3 by a solid 
line. It has indeed a pocket-like molecular shape used by several phenomenological 
approaches [8]. 


3 Matching conditions 


The decay width can be estimated by using the well known expression [8] 


= hu Fint(R) = ~ cons 
T=h = = P(R)y(R) t, (3.1) 


in terms of the Coulomb penetrability and reduced width 


_ _* 
PR) = FR) 
h 
(R) = aun) (3.2) 


2uR ant 
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The general matching condition between the internal and external wave functions 
for a spherical emitter can be written in terms of the logarithmic derivative 


int 2) = Eat) 
Fint(®) Jen A) 


In principle Eqs. (3.1) and (3.3) should not depend upon the radius R, i.e. the 
internal microscopic and external phenomenological solutions should be parallel, at 
least for some region between the touching configuration and Coulomb barrier. This 
common condition, called “plateau condition”, is automatically satisfied when both 
internal and external components are provided as one solution of the Scrédinger 
equation. In our case the internal solution is provided by the amplitude (2.4) and 
the proof of the “plateau condition” is a test of validity for the microscopic approach. 
This task is difficult to achieve for two main reasons. 


(3.3) 


(i) The microscopic estimate of the formation amplitude in terms of two-proton 
and two-neutron orbitals provided by the above described microscopic procedure 
gives a decay width by two orders of magnitude smaller than the experimental 
value. 

(ii) The external logarithmic derivative is proportional to the Coulomb parameter 
x. Indeed, the external solution fe.;() beyond the barrier region is proportional 
to the irregular Coulomb wave Go, because the regular solution practically vanished 
here and a very good estimate (within a 2% error) is given by the WKB ansatz 


Gilp,x) & Golp,x)Ci 


Go(p,x) = (ctg a)? exp [x (a — sin a cos a)| 
Cl = exp oon | ; (3.4) 
x Ve 


where 


2 p Q — QR 


cos a=— 


x Vo(R) 222 


21, (3.5) 


It turns out that the internal formation amplitude, provided by the above described 
microscopic procedure, is insensitive to this parameter. 
4 Hybrid single particle basis 


This drawback can be cured by using a hybrid basis containing an a-cluster compo- 
nent in addition to the standard sp basis Wreljm(X) = (LG ain) x = (r,s), provided 
as an eigenstate of the standard mean field 


Wsp = Wreljm =p ti ehie . (4.1) 
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Figure 4: Two ho potential wells generating the hybrid sp basis (4.1). Solid curve 
generates standard sp orbitals, while the dashed curve is responsible for the cluster 
component [13]. 


This idea was employed in Ref. [14], by using a cluster component of about 30% 
in order to describe the experimental decay width of ?7!*Po. Anyway, the condition 
(ii) was not fulfilled, because only one emitter was analyzed. Notice in this context 
that the microscopic formation amplitude (2.15) contains the ho function, depending 
on four times the ho parameter 3. This is due to the fact that sp proton (7 = p) and 
neutron (7 = 7) orbitals are expanded in terms of a ho basis with the parameter 6 


No 
brajm(®) = > aren [OP () Bx4(8)] (4.2) 
N=2n-+l=0 a 


depending on the radial quantum number n. 
In order to achive both matching conditions (i) and (ii) we expand the cluster 
wave function in a similar way 


Ni 


teus®) = > aen [HP @x1(5)] (4.3) 


gm 
N=2n+l=No+1 


but depending on a different ho parameter {, [{13, 15]. In order to ensure a larger 
radial tail of the sp orbital, providing a larger value of the decay width, the second 
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R (fm) 


Figure 5: Logarithm of the ratio between the theoretical decay width and the exper- 
imental value versus radius if one considers No = 2n+/ = 6 major shells (dot-dashed 
curve), N; = 9 shells with f = 1 (dashed curve) and f = 0.2 (solid curve) for the 
last three major shells [15]. The decay process in 2!? Po +298 Pb + a. 


ho parameter 


Bo=fB , (4.4) 


should be smaller, i.e. f <1, thus corresponding to a larger ho well, as seen in Fig. 
4 (dashed curve). Notice that the solid ho well generates standard sp orbitals. 

In Fig. 5 we analyzed the logarithm of the ratio between the theoretical decay 
width and the experimental value for the decay process in 7!* Po +?°8 Pb+ a versus 
radius by considering No = 2n +1 = 6 major shells in the diagonalisation procedure 
(dot-dashed curve), N; = 9 major shells with f = 1 (dashed curve) and f = 0.2 for 
the last three major shells (solid curve) [15]. In the first two cases the dependence 
versus radius is strong and the decay width underestimated, while in the last case 
the dependence becomes quasiconstant and the experimental width is reproduced in 
the barrier region R € [10,12] fm. 

In Ref. [7] we performed a systematic analysis of a-decay widths between ground 
states of even-even nuclei. We used the pairing approach in order to determine two- 
body X-coefficients. Details are given in Section 5.1. It turns out that more than 
90% of the decay width is given by the cluster part of the expansion (4.1) with 
f ~ 0.7. On the other hand, the condition (ii) is fulfilled if one considers the 
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Figure 6: Proton and neutron potentials in 74°Rn (solid lines) and the corresponding 
Woods-Saxon mean-field potentials (dashed lines) [16]. The cluster parts are as in 
Eq. (4.5) with parameters b, = 1 fm, and r, = 1.3(A!/3 + 41/3) fm. 


relation 6. ~ x, due to the fact that the logarithmic derivative of the internal wave 
function (peaked on the nuclear surface) is proportional to the cluster ho parameter 
Be, while the external wave function is proportional to the Coulomb parameter x. 

Let us stress on the fact that we kept unchanged the parameters of the Woods- 
Saxon potential and changed only the ho parameter for the last shells, describing 
sp orbitals at large distances. This hybrid basis can be generated in a natural way 
by considering the standard Woods-Saxon mean field plus a Gaussian correction on 
the nuclear surface [16] 


V(r) = Wyrs(r) + Vaue(r) 
-(52)| 


ve chavs = Vowel Te) (4.5) 


clus 


lus(T) = Vyexp 


clus 


The potentials for proton and neutrons are plotted in Fig. 6 by solid curves. In Ref. 
[16] it was shown that the Woods-Saxon part of this potential provides the low- 
energy and relatively low-spin states that describe spectroscopic properties, while 
the pocket-like part provides high-energy states with high spin, describing the decay 


process. In Fig. 7 (a) we plotted the proton cluster strength Vis satisfying the 
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Figure 7: (a) Proton clustering strength —V/,, as a function of the neutron number 
for Z < 82, 82 < N < 126 (squares), Z > 82, 82 < N < 126 (triangles), and 
Z > 82, N > 126 (circles) [16]. (b) Logarithm of the reduced width at the touching 
configuration vs neutron number N. (c) The difference between the Q-value and 


Coulomb barrier vs N. 


“plateau condition” 


1 ‘a T(R) 
a ike =i: 4.6 
san | ay MOTH (4.6) 


in the barrier region, versus neutron number for various regions between magic 

numbers. It can be seen that the dependence of the cluster strength is similar to 

the reduced width plotted in the panel (b) and minus the fragmentation potential 
Vprag = —(VB — Q), plotted in the panel (c). 


5 Proton-neutron versus a-like correlations 


In Ref. [12] we analyzed the influence of proton-neutron correlations on a-decay 
in the region above !°°Sn. The spectroscopic factors in this region have indeed the 
largest values. One possible reason is the enhancement of the proton-neutron resid- 
ual interaction, due to the fact that proton and neutron valence orbitals correspond 
to the same major shell. 
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5.1 Formation amplitude 


The amplitude for the decay process of the parent nucleus into one a-cluster and 
a daughter nucleus is given by the overlap integral (2.1). The internal a-particle 
wave function is given by the relative-motion factor of the product between two- 
proton and two-neutron lowest harmonic oscillator states with the parameter Bq ~ 
0.5 fm-? [4]. For sufficiently large inter-fragment distances with a weak overlap of 
nuclear densities, the factorization (2.3) is valid for the parent wave function and the 
formation amplitude simplifies to (2.4). We look for the form of the a-like quartet 
W, in the Fock space representation as 


PL = SO Zpmmn(Lr,L2) [Phr, ® Phaca| 
L1,L2 


a3 S> Zpnpn(L1, La) Phas @ Phira|_ ? (5.1) 
I1,L2 


0 


expressed in terms of the coherent pairs 


Pe estar = we (j1J2; L) PI, (t141 7232) ’ (5.2) 
jij2 
where 
il 
Pe tages = |, all 5.3 
se (m15i72J2) 1+ Orr D1 50 “nh 7 rai LM 


We shall use the fact that the second term of Eq. (5.1), representing the contribution 
of proton-neutron pairs, may be recoupled to a standard proton-proton and neutron- 
neutron form, i.e. 


toma f Spat ee SOTA oe fe at 

Ge S d,,| : leh 8 ot 0 FR »! [leh 8 ty], ® leh, 8 ct, | se AB) 
l 

By using the spherical symmetry, the form of the quartet becomes 


= 1 die: 1 
PL. = Lope > | pO) i let, @d,] > Xmn(I'I's0) Fa ae ad, 
; 7 


l .. ee 
— Zonpn >) x Xpn(I5:0) Xpn(i'9’;0) [et @ a et, @ <, | I, 
jd dam 
(5.5) 


We perform the overlap integral by using the coordinate representation 
let, @et,] + ALdn (x1) © thin (x2)} 


~ = {fj (x1) ® Dio (x2) ear — ar (K2) @ boca} - 
(5.6) 
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We then use the ho representation 


retim(X) = J en(reljm) nh) » (5.7) 
where we considered the sp basis within the jj scheme 
Prim) = [Ou (e) @ xq (9)] 
oD) = RO (”’)iEVim(P) - (5.8) 


By using the recoupling to the LS scheme, the changing from absolute to relative 
and cm proton and neutron pair coodinates, and then by recoupling to relative and 
cm a-coordinates one obtains the following result 


F(R) = 0 W(Na)oN.9(R) 5 (5.9) 
Na 
where 
W(Na) = Zppnn Y> Gp(Np) Gn(Nn)(na0Na0; 0|Np0N 0; 0) ZE8> 
Np,Nn 
Zonpn >.> Gpn(Np; Na, L)(na0No0; 0/NpLNnL;0) LER”) | 
Np Nn 
(5.10) 
The G coefficients contain all recoupling transformations. The ppnn coefficients are 
. 11 Dhol & 
Gr(Nr) = D7 Br(manej) ((U)0(55)05 Ol(l5)5 (U5): 0) 
n1N2Jj 
x S°(-)/(n-0N,0; Ori inal; 0) ree FS 
(5.11) 
while for the pnpn part we have 
Gon(Np» Nn L) = D/L Xpn(G5;0)Xpn(9'H's 0) Byn (95's Np, L) Byn(I5's Nn L) 
Jd 
tLe 1 1 
Wf IA — 0: D(a) BY 
x (GF) LS 5)0: Mig) i i 5)F 5 £) 
(5.12) 


They are expressed in terms of the coefficients 


X77 (9550) Cri (TI) Cn2(7J) 5 
V2 S~ ens (9) Cn (77')tY' (n,ON,L; L[maljnalyys L) TEGO . 


Nr ning 


B, (n1N9; J) 
Bol iF Nex) 


(5.13) 
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5.2 Pairing correlations 


We will consider in our analysis a schematic isovector pairing Hamiltonian with 
constant strengths 


Ce ns 
H= D7 dieiN- DU FDA Phas 
a 


T=p.n j L=pp,nn,pn 
(5.14) 
given in terms of the number of particles and the pair operators of Eq. (5.3) 
_ y y 
Nj — > (CL jmCrjm ir CamCrim) , 
m>0 
P!. = Pl (pjpj) , Pl. = Ply(njnj) , Pl, = Phy(pinj) 


(5.15) 


5.3. pn-BCS treatment 


When taking into account proton-neutron correlations, the quasiparticles carry an 
additional degree of freedom as the full Bogoliubov transformation also mixes pro- 
tons and neutrons 


em tat Ute, 4s PSUS (5.16) 
T=pn 


As a particular case, one recovers the decoupled BCS version for protons k = 1, tT = p 
and neutrons k = 2, rT = n, respectively. The pairing gaps read 


2j+1 
A, = ca 9 Ukr jUkrj > TH=pPi,n, 
kj 
29 ae Uknj Vkpj + UkpjVknj 
kj 
(5:17) 
The coherent BCS pair amplitudes of Eq. (5.2) are given by 
1 V2ag4+1 
Xj = xP ler,j¢r,Jlp |P) = a S WieniPirj > THPin, 
k 
1 2j3+1 UknjVkpj + UkpjVkv4 
Xpnj = 3D [epgenglg |P) = ers eg : 
k 
(5.18) 
For the a-like quartet to be isoscalar we require 
2 1 
Zopnn — V3 ) Zynpn = VB : (5.19) 
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5.4 Exact solutions 


In the case of single species pairing (Gpn = 0), we may easily find the form of the 
ground state by diagonalising the Hamiltonian in the 2n-particle basis. A generic 
state is written as 


is Sa a PO (5.20) 


iy>...>tn 


t_ tte ot 
where P! = a [el © all 
As far as the static properties are concerned, the BCS and diagonalisation treat- 


ments are shown below to agree well. However, the cohererent pair amplitudes X;,; 
are generally overestimated by the BCS treatment. 


In Fig. 8 we plotted the exact diagonalization results for Z = N = 52 by using a 
decoupled basis for proton and neutrons (dashes) and coupled pn basis (solid line). 
These results are indeed larger in comparison with the BCS (dashes) and pn-BCS 
cases (thick solid line), but still they underestimate the experimental decay widths. 
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& 
= 10-1074 
“" 
4 
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Figure 8: Comparison of formation amplitudes in various approaches. The proton 
spectrum is considered identical to the neutron spectrum [12]. 


As we already pointed out, the standard mean-field approach cannot supply 
enough four-body correlations to form an a-particle on the nuclear surface, because 
the computed decay widths are too small. One convenient way to simulate the 
missing four-body correlations was described in the previous Section by considering 
an additional Gaussian pocket-like attractive component on the surface region of 
the single-particle mean field [16]. We considered such a component in both proton 
and neutron mean fields with the same parameters 
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Figure 9: The strength of the a-clustering additional part of the mean field potential 
Vy versus the reduced width [12]. 


By adjusting the strength of this component V, , it is possible to reproduce the 
experimental decay widths. In Fig. 9, we plotted this strength within the BCS ap- 
proach versus the experimental reduced width for the region above !°°Sn (triangles) 
and above ?°8Pb (solid circles). One remarks that these values are proportional, 
due to the fact that they were simultaneously determined. Thus, the strength of 
the additional potential is greater for lighter nuclei with Z ~ N. Let us mention 
that we obtained similar values for light nuclei by using the pn-BCS approach (open 
squares). Therefore the pn correlations are not able to explain larger reduced widths 
for N ~ Z nuclei above !°°Sn. The fitted values plotted by the corresponding lines in 
Fig. 9 have similar slopes and intercepts. This feature proves the universal character 
of the a-clustering across the nuclear chart. 


6 Conclusions 


We explained the a-like rotational bands in the 4°Ca+a system by computing the 
a-formation amplitude in terms of single particle Gamow resonances. We demon- 
strated that the systematics of decay widths between ground states cannot be ex- 
plained in terms of a collective four-body state built from two-proton and two- 
neutron orbitals of the standard nuclear mean field. It can be understood only in 
terms of a preformed a-cluster existing on the nuclear surface in addition to the 
above mentioned mean field cluster. Therefore, we used a semi-microscopic hybrid 
model combining the mean field formation with a preformed a-cluster in order to 
explain the order of magnitude of the experimental decay width. We compared 
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proton-neutron with a-like correlations above the double magic nucleus !°°Sn. Fi- 
naly we evidenced the universal character of a-like correlations. 
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